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A widely adopted approach to solving constraint satisfaction problems combines systematic tree 
search with various degrees of constraint propagation for pruning the search space. One common 
technique to improve the execution efficiency is to add redundant constraints, which are con- 
straints logically implied by others in the problem model. However, some redundant constraints 
are propagation redundant and hence do not contribute additional propagation information to 
the constraint solver. Redundant constraints arise naturally in the process of redundant mod- 
eling where two models of the same problem are connected and combined through channeling 
constraints. In this paper, we give general theorems for proving propagation redundancy of one 
constraint with respect to channeling constraints and constraints in the other model. We illustrate, 
on problems from CSPlib (http://www.csplib.org/), how detecting and removing propagation 
redundant constraints in redundant modeling can significantly speed up constraint solving. 
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and Features — Constraints; F.4.1 [Mathematical Logic and Formal Languages]: Mathemat- 
ical Logic — Logic and constraint programming 

General Terms: Performance, Theory 

Additional Key Words and Phrases: Constraint Propagation, Redundant Constraints, Redundant 
Modeling 



1. INTRODUCTION 

Finite domain constraint programming combines backtracking tree search with con- 
straint propagation to solve constraint satisfaction problems (CSPs) [Mackworth 
1977]. Constraint propagation removes infeasible values from the domains of vari- 
ables to reduce the search space. This propagation-based constraint solving frame- 
work is realized in modern constraint programming systems such as ECL^PS*^ [Chea- 
dle et al. 2003], ILOG Solver [1999], and SICStus Prolog [2003], which have been 
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successfully applied to many real-life industrial applications. 

There is usually more than one way of modeling a problem as a CSP. By mod- 
cling a problem as a CSP, we mean the process of determining the variables, the 
associated domains of the variables, and the expressions of the constraints. Find- 
ing a good model of a CSP is a challenging task. A modeler must specify a set of 
constraints that capture the definitions of the problem, but this is not enough. The 
model should also have strong propagation: that is, it should be able to quickly 
reduce the domains of the variables of the problem. Moreover, the implementation 
of propagators to perform constraint propagation should be efficient. Last but not 
least, the model should have a small search space. ^ 

A common technique to increase propagation strength is to add redundant con- 
straints, which arc logically implied by the constraints of the model. Adding re- 
dundant constraints can be beneficial since the constraint solver may extract more 
information from these redundant constraints. However, some logically redundant 
constraints are propagation redundant, and hence do not contribute additional prop- 
agation information to the constraint solver. Generally, we only want to add redun- 
dant constraints that are propagation non-redundant to reduce the search space. 

An important source of logically redundant constraints is in redundant model- 
ing [Cheng et al. 1999]. A problem can be modeled differently from two viewpoints 
using two different sets of variables. By connecting the two different models with 
channeling constraints, which relate valuations in the two different models, stronger 
propagation behavior can be achieved in the combined model. However, the addi- 
tional variables and constraints impose extra computation overhead. Given each 
model is complete and only admits the solutions of the problem then each model is 
logically redundant with respect to the other model plus the channeling constraints. 
In many cases, some of the constraints are also propagation redundant with respect 
to the other constraints in the combined model. By reasoning about propagation 
redundancy, we can improve redundant modeling by just keeping the constraints 
which give beneficial new propagation. 

In this paper, we introduce the notion of restrictive and unrestrictive channel 
functions to characterize channeling constraints. We study the propagation be- 
havior of constraints based on the notion of propagation rules, which capture each 
possible propagation by a constraint. This allows us to systematically determine if 
a propagator is redundant with respect to the propagators of a set of constraints 
through simple implication tests. We give general theorems for proving propagation 
redundancy of constraints involved in redundant models. We focus on propagators 
that perform (the combination of) two popular propagation techniques, namely do- 
main propagation [Van Hentenryck ct al. 1998] and set bounds propagation [Gervet 
1997]. We illustrate, on problems from CSPLib (http://www.csplib.org/), how 
detecting and removing propagation redundant constraints can significantly speed 
up solving behavior. This paper is a revised and extended version of our earlier 
work [Choi ct al. 2003a; 2003b]. 

The remainder of the paper is organized as follows. In Section 2 we introduce 
propagation-based constraint solving and propagation rules, a way of enumerating 
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the different propagation behaviors of a propagator. In Section 3 we define a broad 
form of channeling constraints that are covered by our approach. In Section 4 we 
give theorems that allow us to show which constraints in a redundant model are 
not causing extra propagation and can be removed. In Section 5 we discuss some 
important issues in studying about propagation redundancy. In Section 6 we give 
experimental results showing the benefits of detecting and removing propagation 
redundant constraints. In Section 7 we discuss related work. In Section 8 we 
summarize our contributions and shed light on future directions of research. 

2. BACKGROUND 

In this paper we consider integer and set constraint solving with constraint propa- 
gation and tree search. Boolean constraint solving is considered as a special case of 
integer constraint solving. Our notations, although different from the conventional 
CSP literatures, allow us to express the theoretical framework in a simpler manner. 

2.1 Variables and Domains 

We consider a typed set of variables V = Vj U Vs made up of integer variables Vi, 

for which we use lower case letters such as x and y, and sets of integers variables 
Vs, for which we use upper case letters such as S and T. We use v to denote 
variables of either kind. 

Each variable is associated with a finite set of possible values, defined by the 
domain of the CSP. A domain D is a complete mapping from a fixed (countable) 
set of variables V to finite sets of integers (for the integer variables in V/) and to 
finite sets of finite sets of integers (for the set variables in Vs)- A false domain D is 
a domain with D{v) — for some v. A singleton domain D is such that |-D(t))| = 1 
for all V € V. The intersection of two domains Di and D2, denoted Di n D2, is 
defined by the domain D^lv) = Di{v) Cl D2{v) for all v. A domain Di is stronger 
than a domain D2, written Di C D2, if Di{v) C D2{v) for all variables v. A 
domain Di is equal to a domain D2, denoted D\ = D2, if Di{v) = D2{v) for all 
variables v. We also use range notation whenever possible: [/..u] denotes the set 
{d \ I < d < u} when I and u are integers, while [L ..U] denotes the set of sets of 
integers {A \ L C A C U} when L and U are sets of integers. We shall be interested 
in the notion of an initial domain, which we denote Dinit ■ The initial domain gives 
the initial values possible for each variable. In effect an initial domain allows us to 
restrict attention to domains D such that D C Dinu . 

2.2 Valuations, Infima and Suprema 

A valuation is a mapping of integer variables {xi € V/) to integer values and 
set variables {Si G Vs) to sets of integer values, written {xi ^ di,...,Xn 1— > 
dn, Si H- > Ai, . . . , Sm Am} where di G D{xi) and Aj G D{Sj). Let vars be the 
function that returns the set of variables appearing in an expression, constraint or 
valuation. Given an expression e, 6{e) is obtained by replacing each v G vars{e) 
by 9{v) and calculating the value of the resulting variable free expression. In an 
abuse of notation, we define a valuation 9 to be an element of a domain D, written 
6* G D, if 9{vi) G D{vi) for aU Vi G vars{9). 

Define the infimum and supremum of an expression e with respect to a domain 
D as mine = inf{^(e)|^ G D} and sup£,e = sup{^(e)|^ G D}. The ordering ^ 
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used by inf and sup depends on the type of the expression. If e has integer type 
then di ^ d2 iff di < d2, while if e has set of integer type then di ^ ^2 iff di C ^2. 
Note that these values may not exist for arbitrary domains and set of integer type 
expressions. Later we shall restrict ourselves to domains and expression where 
infimum and supremum always do exist. 

2.3 Constraints 

A constraint places restriction on the allowable values for a set of variables and 
is usually written in well understood mathematical syntax. More formally, a con- 
straint c is a relation expressed using available function and relation symbols in a 
specific constraint language. For the purpose of this paper, we assume the usual 
integer interpretation of arithmetic constraints, set operators such as G and C, and 
logical operators such as -i, A, V, and A CSP consists of a set of constraints 
read as conjunction. We use true to denote the empty (always satisfiable) conjunc- 
tion of constraints. We define solns{c) = {6 \ vars{6) = vars{c) A \=e c}, that is 
the set of 9 that make the constraint c hold true. We call solns{c) the solutions 
of c. In some cases, constraints can also be defined directly by giving the set (or 
table) solns{c). We sometimes treat an integer constraint c as an expression with 
value 1 if true and if false. We can understand a domain D as a constraint in the 
obvious way, D Avev VdeD(w) v = d. A constraint c is logically redundant 

with respect to a constraint c' if |= c' — > c, that is c holds whenever c' holds. 

2.4 Propagators and Propagation Solvers 

In the context of propagation-based constraint solving, a constraint specifies a prop- 
agator, which gives the basic units of propagation. A propagator / is a monoton- 
ically decreasing function from domains to domains, i.e. Di C D2 implies that 
f{Di) E /(-D2), and f{D) C D. A propagator / is correct for constraint c iff for 
all domains D 

{0\0GD}n solns{c) = {e\0G f{D)} n solns{c) 

This is a weak restriction since for example, the identity propagator is correct for 
all constraints c. We will usually assume that a propagator / for constraint c is 
checking, that is, if Z) is a singleton domain, then f{D) = D iS30 £ Df] solns{c). 
A checking propagator correctly determines the satisfiability of the constraint c for 
singleton domains. 

A propagation solver for a set of propagators F and current domain D, solv{F, D), 
repeatedly applies all the propagators in F starting from domain D until there is no 
further change in resulting domain. solv{F, D) is the largest domain D' Q D which 
is a fixpoint (i.e. /(£)') = D') for all f € F. In other words, solv{F,D) returns a 
new domain defined by 

iter{F,D) = \1 f{D) and solv{F,D) = gi^{\d.iter{F,d)){D) 

where gfp denotes the greatest fixpoint w.r.t C lifted to functions. 
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2.5 Domain and Set Bounds Propagators 

Propagators are often (but not always) linked to implementing some notion of local 
consistency. The most well-studies consistency notion is arc consistency [Mack- 
worth 1977] which ensures that for each constraint, every value in the domain of 
the first variable, has a supporting value in the domain of the second variable which 
satisfied the constraint. Arc consistency can be naturally extended to constraints 
of more than two variables. This extension has been called generalized arc consis- 
tency [Mohr and Masini 1998] , as well as domain consistency [Van Hentenryck et al. 
1998] (which is the terminology we will use), and hyper-arc consistency [Marriott 
and Stuckey 1998]. 

A domain D is domain consistent for a constraint c if D is the least domain 
containing all solutions 9 G D oi c, i.e, there does not exist D' \Z D such that 
6 G D A 6 G solns{c) 6 £ D'. A set of propagators F implement domain 
consistency for a constraint c, if solv{F, D) is always domain consistent for c. Define 
the domain propagator for a constraint c as 



Example 2.1. Consider the constraint c = xi = 3x2 + 5.r3. Suppose do- 
main D{xi) = {2,3,4,5,6,7}, D{x2) = {0,1,2}, and ^(xs) = {-1,0,1,2}. The 
solutions ^ G D of c are 6i = {xi i-^ 3,X2 ^,X3 i-^ 0}, 02 = {xi ^ 5,X2 
0,0:3 1-^- 1}, and 6*3 = {xi 1-^ 6,X2 2,2:3 0}. Hence, dom{c){D) — D' where 
D'{xi) = {3, 5, 6}, D'{x2) = {0, 1, 2}, and D'{xz) = {0, 1}. D' is domain consistent 
with respect to c. □ 

Set bounds propagation [Gervet 1997] is typically used where a domain maps a 

set variable to a lower bound set of integers and an upper bound set of integers. 
We shall enforce this by restricting our attention to domains where the D{S) is 
a range, that is D{S) = {A \ infD(S') C ^ C supj5(S')}. This is managed by 
only using set bounds propagators, which maintain this property. The set bounds 
propagator returns the smallest set range which includes the result returned by 
the domain propagator. Define the set hounds propagator for a constraint c where 
vars{c) C Vs as 



Example 2.2. Consider the constraint c = 5*1 C S2, suppose the domain 



D{Si) = [{1}.. {1,2, 3, 4}], D{S2) = [0..{1,2,3}]. Then, D' = sb{c){D) where 
£>'(5i) = i?'(52) = [{l}..{l,2,3}]. □ 



A constraint can involve both integer and set variables. We define the domain 
and set bounds propagator dsb{c) for a constraint c as follows: 



Note that as defined dsb{c) = dom{c) when vars{c) C V/. Prom now on we shall 
restrict attention to dsh propagators. 
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where v S vars{c) 
otherwise 




where v £ vars{c) fl Vs 
otherwise 



Example 2.3. Consider the domain propagator dom (151 — x), applied to a 
domain D where D{x) = {2}, and D{S) = [0 .. {1, 5, 8} ]. Then the solutions 
9 e D are {x ^ 2, S ^ {1,5}}, {x ^ 2,S ^ {1,8}}, {x ^ 2,S ^ {5,8}}. So 
D' = dom{\S\ = x){D) gives D'{S) = {{1, 5}, {1, 8}, {5, 8}}. The domain and set 
bounds propagator rfs6(|S'| = x) instead determines D = dsb{\S\ = x){D) since 
n{{l, 5}, {1, 8}, {5, 8}} = and U{{1, 5}, {1, 8}, {5, 8}} = {1, 5, 8}. □ 

2.6 Atomic Constraints and Propagation Rules 

An atomic constraint represents the basic changes in domain that occur during 
propagation: the elimination of a value from an integer domain, or the addition of 
a value to a lower bound, or removal of a value from an upper bound. An atomic 
constraint is one of Xi = d, Xi ^ d, d G Si or d ^ Si where Xi gVi, d is an integer, 
and Si e Vs. 2 

A propagation rule is of the form C ^ c where C is a conjunction of atomic 
constraints, c is an atomic constraint, and ^ C — > c. For notational convenience we 
shall write extended rules C ^ C where C is a conjunction of atomic constraints 
as a shorthand for the set of rules {C ^ c \ c G C'}. 

A propagation rule C ^ c defines a propagator (for which we use the same 
notation) in the obvious way. 

N/r^N/ N ( \6(v) \ 9 & D A9 e solns(c)} if mrs(c) = {v} and D ^ C 

Wc can characterize an arbitrary propagator / in terms of the propagation rules 
that it implements. A propagator / implements a propagation rule C >^ c if for 
each D C Dinit whenever |= D — > C, then 1= f{D) — » c. Let rules{f ) be the set of 
rules implemented by /. Tliis definition of / is often unreasonably large. In order 
to reason more effectively about propagation rules for a given propagator, we need 
to have a minimal representation.^ Let prop{f) C rules{f) be a set of propagation 
rules such that solv(prop{f), D) = solv{rules{f), D) for all D C Dinit, and there 
does not exist a prop{f) Cprop{f) for which solv{prop{f), D) = solv{rules{f), D) 
for all D C Dinit- That is, all propagation caused by / is also caused by the 
minimal set of propagation rules prop{f). Notice that prop{f) is not unique. 

Example 2.4. For Boolean constraint c = 212 + ^13 + 214 — 1 where Dinit{zi2) = 
Dinitizis) = Drmiizii) = [0.. 1] thc propagation rules prop{dsh{c)) arc 

212 = 1 ^ 213 = 0, 214 = 212 = 0, 213 = ^ 2i4 = 1 

213 = 1 ^ -212 = 0, 214 = 212 = 0, 2i4 = ^ 2i3 = 1 

214 = 1 ^ 212 = 0, 213 = 213 = 0, 2i4 = ^ 2i2 = 1 

Note that propagation rules prop{f) for propagators / for Boolean constraints need 
only ever involve equations since |= Dinit {z — b ^ z ^ I — b) for b G {0, 1}. □ 

A key lemma for domain and set bounds propagators dsb{c'), is that the propa- 
gation rules implemented are exactly those C ^ c where c' implies C ^ c. 



^Atomic constraints of the form Xi = d are not strictly necessary for propagation rules. They are 

equivalent to removing all other values from thc domain. 

^Both Brand [2003] and Abdennadher and Rigotti [2002] give effective methods for creating min- 
imal representations of any constraints in terms of propagation rules. 
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Lemma 2.5. Given a constraint c' , dsb{c') implements C ^ c iff \= {Dinu A 
c')-(C-c). 

Proof. (=>) Suppose dsh{c') implements C ^ c. Then for each D C Dinu we 
have that \= D ^ C means that |= dsb{c'){D) — > c. 

Suppose to the contrary ^ {Dinu Ac') (C c). Thus there exists solution 
6 G Dinit that satisfies c' A C A -■(c). We build a domain Dg from a valuation 6 
as follows. Dg{v) = {0{v)} for all u G vars{6), while Dg{v) = Dinniv) otherwise. 
Now \= Dg ^ C and dsb{c'){Dg) = Dg since 6* is a solution of c'. But then 
^ dsb{c'){Dg) c. Contradiction. 

(<;=) Suppose 1= {Dinit A c') — ^ (C ^ c). Then for any domain D C -Dimt where 
^ D ^ C, we have that \= {D A c') ^ c since |= D — > Dmu- Hence, every solution 
^ e of c' is a solution of c. 

Suppose to the contrary dsb{c') docs not implement C ^ c. Then there exists 
D E D,mt, such that ^ D ^ C and D' = dsb{c'){D) but ^ £)' c. We now prove 
a contradiction for each form of c: 

— If c = d G S, then ^ D' — > c means d ^ infDi{S). Hence by the definition of 
dsb(c'), there is a solution 6 G D oi c' where d ^ 0{S). Contradiction since 6 is 
not a solution of c. 

— If c = d ^ S, then ^ D' — > c means d e sup£)/(S'). Hence by the definition of 
dsb(c') there is a solution 6 G D of c' where d € 0{S). Contradiction since 6 is 
not a solution of c. 

— If c = X ^ d, then ^ £)' — > c means d e D'{x). Hence by the definition of dsb{c'), 
there is a solution 6 G D of cs where 6{x) = d. Contradiction since 6 is not a 
solution of c. 

— If c = X — d, then ^ D' ^ c means d' £ D'{x) where d' ^ d. Hence by the 
definition of dsb{c'), there is a solution 6 G D oid where 6{x) = d' . Contradiction 
since 9 is not a solution of c. 

□ 

3. CHANNELING CONSTRAINTS 

Redundant modeling [Cheng et al. 1999] models a problem from more than one 
viewpoint. By joining two models using channeling constraints, we can get the 
advantage of both sources of propagation. 

Assume we have one model of the problem Mx using variables X, and another 
model My using disjoint variables Y. Channeling constraints can be used to join 
these two models together by relating X and Y. There is no real agreement, as yet, 
as to precisely what channeling constraints are. For the purposes of our theorems 
we define a channeling constraint as follows. 

Let Ax be the atomic constraints for Dinit on variables X, and Ay be the atomic 
constraints for D„iu on variables Y. A channel function 'O' is a bijection from atomic 
constraints Ax to Ay- 

A channeling constraint (or simply channel) Co is the constraint 

A (c^ow) 

c^Ax 
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The channel propagator is the set of propagation rules inferred from the 
channel function <(). 

^0 = U {c ^ 0(c), 0(c) ^c} 

ceAx 

Note for channel function 0, by definition 0^^ is also a channel function, and 
and C<^-i, as well as F<^ and F^-i, are identical. 

We now illustrate how common channels fit into this framework. 

3.1 Permutation Channels 

A common form of redundant modeling is when we consider two viewpoints to 
a perrnAitation problem [Geelen 1992]. In a permutation problem, the objective 
is to find a bipartite matching between two sets of objects A = {oi, . . . , a„} and 
B = {bi, . . . , bn} satisfying all other problem specific constraints. Generally, we can 
model a permutation problem from two different viewpoints. In the first viewpoint, 
we assign objects from B to A. We use the set of variables X = {xi, . . . ,a;„} to 
denote objects in A, and the domain D{xi) = {1, . . . , n}, for all 1 < i < n to denote 
objects in B. The second viewpoint swaps the role between A and B, i.e. assign 
objects from A to B. We use the set of variables Y ~ {yi, ...,?/„} to denote objects 
in B, and the domain D{yj) = {1, . . . , n}, for all 1 < j < n to denote objects in A. 

The permutation channel function 1X1 is defined as N (.t; = j) = (yj = i) and 
^{xi 7^ j) = ivj 7^ i) for all 1 < i,j < n. The permutation channel Cm is equivalent 
to the conjunction of constraints 

n n 

/\ /\{xi=j Vj = i) 
i=i j=i 

Example 3.1. Langford's Problem The problem "prob024" in CSPLib is an 

example of permutation problem. The problem is to find an (m x n)-digit sequence 
which includes the digits 1 to n, with each digit occurring m times. There is one 
digit between any consecutive pair of the digit 1, two digits between any consecutive 
pair of the digit 2, . . . , and n digits between any consecutive pair of the digit n. 

Smith [2000] suggests two ways to model the Langford's problem. We use the 
(3 X 9) instance to illustrate the two models. In the first model, Mx, we use 27 
variables X — {a;i, . . . , X27}^ which we can think of as li, I2, 13, 2i, . . . , 92, 93. Here, 
li represents the first digit 1 in the sequence, I2 represents the second digit 1, and 
so on. The initial domain of these variables, Dinit{xi) = {1, . . . , 27} for 1 < i < 27, 
represents the positions of the digit Xi in the sequence. We enlist the constraints 
of Smith's model as follows: 

— (LXl) disequality constraints: VI < i < j < 27. Xi ^ Xj 

— (LX2.1) separation constraints: VI < i < 9. Xsi-i = Xsi-2 + + 

— (LX2.2) separation constraints: VI < i < 9. Xsi = x^i-i + + 

In the second model, My, we again use 27 variables Y = {yi, . . . , 2/27} to represent 
each position in the sequence. The initial domain of these variables, Dinit{yi) = 
{1,...,27} for 1 < i < 27, corresponds to the digits li, I2, 13, 2i, . . . , 92, 93 in 
position yi of the sequence. The constraints are: 

8 



— (LYl) disequality constraints: VI < i < j < 27. yi ^ yj 

— (LY2.1) separation constraints: VI < i < 9.V1 < j < 27 - 2{i + 1). = 3i - 2 <^ 
= 3i - 1 

— (LY2.2) separation constraints: VI < i < 9.V1 <j<27- 2{i + 1). yj = 3i-2-^ 

yj+2{i+i) = 3j 

— (LY3) separation constraints: VI < i < 9.V(28 - 2{i + 1)) <j< 27. yj + Zi-2 

The permutation channel for the two models is simply Xi = j <^ yj = i for all 
l<i,j<27. □ 

Example 3.2. All Interval Series Problem The problem "prob007" in CSPLib 
is from musical composition. The problem is to find a permutation of n numbers 
from 1 to n, such that the differences between adjacent numbers form a permutation 
from 1 to n — 1. We give two ways to model the problem. The first model slightly 
modifies the model suggested by Puget and Regin [2001], and the the second model 
slightly modifies the model suggested by Choi and Lcc [2002]. 

The first model, Mx, consists of n variables, X = {xi, . . . , .t„}. Each Xi denotes 
the number in position i, and Dinit{xi) = [ 1 .. n] for 1 < i < n. We introduce auxil- 
iary variables, {wi, . . . , Un-i} that denote the difference between adjacent numbers, 
where Dinit{ui) = [ 1 .. n — 1 ] for 1 < i < n — 1. The constraints are: 

— (IXl.l) disequality constraints: Ml <i < j <n. Xi ^ Xj 
— (1X1.2) disequality constraints: VI < i < j < n — 1. Ui ^ Uj 
— (1X2) interval constraints: VI < i < n — 1. Ui = \xi — Xi+i\ 

The second model. My, also consists of n variables, Y = {t/i, . . . ,y„}. Each y^ 
denotes the position for the number i, and Dinit{yi) = [l..n]forl<i<n. The 
auxiliary variables {wi, . . . , Vn-i} denote the position where the difference value of 
1 to n — 1 belongs, and Dinit (vi) = [ 1 .. n — 1 ] for 1 < i < n — 1. The constraints 

are: 

— (lYl.l) disequality constraints: VI < i < j < n. yi ^ yj 

— (IY1.2) disequality constraints: VI < i < j < n — 1. Vi Vj 

— (IY2.1) interval constraints: VI < j < j < n. (j/j — yj = 1) => {vj-i = yj) 

— (IY2.2) interval constraints: VI < i < j < n. {yj — t/j = 1) {vj-i = t/j) 

The (IY2.1) and (IY2.2) constraints enforce that if yi and yj are adjacent, the 
position for their difference must be the smaller of them. In the second model, 
observe the fact that only the numbers 1 and n can give us the difference of n — 1. 
Therefore, we can add the following redundant constraints: 

(IY3): (|yi - y„| = 1) A (w„_i = min(t/i, y„)), 

which requires yi and y„ to be adjacent. 

The permutation channels for this problem arc more interesting because we have 
two distinct kinds of variables in each model, each of which is related by a per- 
mutation channel. The channels are Xi = j 'i?^ yj = i for all 1 < < n and 
Ui = j <^ Vj = i for all 1 < i, j < n — 1. □ 
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3.2 Boolean Channels 

Another common form of redundant modeling is when we give both an integer 
and Boolean models. Suppose we have an integer model using the integer variables 
X = {xi, . . . ,Xn} and the domain {xi) = [1 ..k]. We can have a corresponding 
Boolean model using the Boolean variables Z = {zij \ 1 < i < n,l < j < k}. Each 
variable encodes the proposition that Xi = j. 

The Boolean channel function A is defined as A{xi = j) = {zij = 1) and A{xi ^ 
j) = {zij — 0) for all 1 < i < n, 1 < j < fc. Note that the atomic constraints 
Zij 7^ 1 and Zij ^ are not needed for Boolean variables since they arc equivalent 
(respectively) to Zij = and Zij = 1. The Boolean channel Ca is equivalent to the 
conjunction of constraints 

n k 

f\ /\iXi=j^ Zij = 1) 

i=lj = l 

Example 3.3. n-Queens Problem The problem is to place n queens on an 
n X n chess board so that no two queens can attack each other. There are two 
common ways to model this problem, i.e., an integer model and a Boolean model. 

The integer model, Mx, consists of n variables, X = {xi, . . . ,Xn}- Each Xi 
denotes the column position of the queen on row i, and Dinit{xi) = {!,... ,n}, for 
1 <i <n. The constraints are: 

— (QXl) column constraints: VI < i < j < n. Xi ^ Xj 

— (QX2.1) diagonal constraints: VI < i < j < n. Xi — i ^ Xj — j 

— (QX2.2) diagonal constraints: VI < i < j < n. Xi + i xj + j 

The Boolean model, Mz, consists of n x n Boolean variables, Z = {zn, . . . , zin, 
. . . , Zni, ■ ■ ■ , Znn}- Each Boolcau variable Zij denotes whether we have a queen at 
row i column j or not. The constraints are: 

— (QZl) row constraints: VI < i < n. ^ij = 1 

— (QZ2) column constraints: VI < j < n. J2^=i ^ij = 1 

— (QZ3.1) main diagonal constraint: X^^^ Za < 1 

— (QZ3.2) main diagonal constraint: Zi{n-i+i) < 1 

— (QZ4.1) diagonal constraints: VI < A; < n — 1. ^j{j+k) < 1 

— (QZ4.2) diagonal constraints: VI < A; < n — 1. Z(^jj^k)j ^ 1 

— (QZ4.3) diagonal constraints: VI < A; < n — 1. Zj(n-j-k+i) < 1 

— (QZ4.4) diagonal constraints: VI < A; < n — 1. z^^k-^^n-j+i) < 1 

We combine the two models using the Boolean channel Xi = j ^ Zij = 1 for all 
l<i<n,l<j <k. □ 

3.3 Set Channels 

Another common form of redundant modeling is where one model deals with integer 
variables, and the other with variables over finite sets of integers, and the relation 
Xi = j holds iff i e Sj. This generalizes the permutation problem to where two or 
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more integer variables can take the same value. Suppose the integer variables are 
X = {xi, . . . , Xn}, where Dinit{xi) = [ 1 .. fc] for all 1 < i < n, and the set variables 
are S = {Si, ...,Sk} where D„„f (S'j) = [0 .. {1, . . . , n} ] for all 1 < j < fc. 

The set channel function {} is defined as {}{xi = j) = {i & Sj) and {}{xi ^ j) = 
{i ^ Sj) for all 1 < i < n, 1 < j < fc. The set channel C{} is equivalent to 

n k 

Example 3.4. Social Golfers Problem The problem "probOlO" in CSPLib is 

to arrange n ~ g x s players into g groups of s players each week, playing for w 
weeks, so that no two players play in the same group twice. Smith [2001] suggests 
two ways to model this problem. 

In the first model we use variables X = {a;/fe|l < Z < n, 1 < fc < w} to denote 
the group which player I plays on week fc, and Dinit{xik) = [1 ..5] for all 1 < Z < 
n,l < k < w. 

The constraints of the problem are expressed as: 

— (GXl) each group has s players: VI < i < g.Vl < fc < w. T,f^-^^{x^k = i) = s 
— (GX2) two players only play in the same group in one week: 

VI < fci ^ fc2 < w.yi <hi=l2<n. -'{xi^ki = xi^ki A xi^k2 = a^/afe) 

The second model uses set variables S = {Sik\l < * < 5, 1 < fc < to denote 
the set of players play in group i on week fc. and Dinit{Sik) = [0 .. {1, ■ ■ ■ in]] for 
all 1 < i < 5, 1 < fc < w. The constraints are expressed as: 

— (GSl) no groups in the same week have a player in common: 

VI < fc < w.Vl <hi^i2<9. Si,k n Si^k = 

— (GS2) each group has s players: VI < i < g.Vl <k <w. \Sik\ = s 
— (GS3) no different groups have more than one player in common: 

VI < zi ^ i2 < 5-Vl < fci ^ fc2 < \Si,k, n Si,k, I < 1 

We can use the set channels to combine the two models, xik = i I & Sik for all 
I <l <n,l < k < w, 1 <i<g. □ 

Example 3.5. Balanced Academic Curriculum Problem The problem, 
listed as "prob030" in CSPLib, is to design an academic curriculum aiming to 

balance the loads in each academic period. Following the description in Hnich et 
al. [2002], we can have both an integer model Mx and set model Ms- 

Given m courses, and n periods, a, b are the minimum and maximum academic 
load allowed per period, c, d are the minimum and maximum number of courses 
allowed per period, ti specifics the number of credits for course i, and J? is a set of 
prerequisite pairs specifying that course i must be taken before course j. 

We introduce a set of auxiliary variables Ij, which is shared by both models, 
to represent the academic load in period j as well as a variable u representing 
the maximum academic load in any period, i.e. u = max{lj | 1 < j < n}. The 
objective function simply minimizes u. We also introduce another set of shared 
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auxiliary variables qj to represent the number of courses assigned to a period. We 
have Dinu{u) = Dinu{lj) = [0 .. S^^ij ] and Dinit{qj) = [ 1 •• m]. 
We have the following constraints that are common to both models: 

— (Bl.l) load allowed per period: VI < j < n. a < Ij < b 

— (B1.2) number of courses allowed per period: VI < j < n. c < Qj < d 

We also add the following redundant constraints: 

— (B2.1) all the credits must be fulfilled: iY.]Li h) = (ElLi^i) 
— (B2.2) all the courses must be taken: (Y^^j^i Ij) ~ ^ 

In the integer model, Mx, the variables X = {xi\l < i < m} represent the 
period to which course i is assigned and Dinit{xi) = [ 1 .. n] for all 1 < i < m. The 
constraints for the integer model Mx are: 

— (BXl) Ij is the load taken in period j: VI < j < n. (X^™ i((a;i = j) x ti)) = Ij 
— (BX2) Qj is the number of courses in period j: VI < j < n. {^^i{xi = j)) = qj 
— (BX3) courses are taken respecting prerequisites: (z R. Xi < Xj 

In the set model, the set variables S = {Sj\l < j < n} represent the set of 
courses assigned to period j and Dinu{Sj) = [ .. {1, . . . , m} ] for all 1 < i < n. 
The constraints for the set model Ms are: 

— (BSl) No course is taken twice: VI < i < j < n. SiH Sj = (l> 

— (BS2) Ij is the load in period j: VI < j < n. {J^ieS- ~ 

— (BS3) qj is the number or courses in period j: VI < j < n. \Sj\ = qj 

— (BS4) courses are taken respecting prerequisites: 

e R.yi <k<n- l.Vl <k' <k. {i G Sk) =^ {j ^ Sk') 

We can use the set channels to combine the two models, xt = j <^ i & Sj for all 
l<i<m,l<j <n □ 

3.4 Channels between Set and Boolean Models 

A very uncommon form of redundant modeling is when wc give a set model and 
a Boolean version of this model. The reason it is uncommon is that there is no 
natural gain in expressiveness in moving to the Boolean model. 

Suppose the set variables arc {Si, . . . , iSfc}. where D„iit{Si) = [0 .. {1, . . . , n} ], 
and the Boolean variables are Zij,l < i < k,l < j < n. The set2bool channel 
function — is defined as —{j S Si) = {zij = 1) and —{j ^ Si) = {zij = 0). The 
set2bool channel is equivalent to 

k n 

/\/\(je5,^Zy = l) 

With — channel, wc can map common set constraints (c) to Boolean constraints 
(—(c)) as given in Figure 1. We shall prove that set bounds propagation of set 
constraints (c) is equivalent to domain propagation for the corresponding Boolean 
constraints (—(c)). 
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Fig. 1. Mapping of Common Set Constraints to Boolean Constraints 
4. PROPAGATION REDUNDANT CONSTRAINTS 

We shall be interested in reasoning about redundancy with respect to sets of prop- 
agators. We say a set of propagators Fi is stronger than a set of propagators F2-, 
written Fx '> F2, if solv{Fi,D) C solv{F2,D) for all domains D C We say 

a set of propagators Fi is equivalent to a set of propagators F2, written Fi ~ F2, 
if solv{Fx,D) = solv{F2,D) for all domains D C Dinit- A propagator / is made 
propagation redundant by a set of propagators F if F » {/}. Our main aim is to 
discover and eliminate propagation redundant constraints and/or propagators. 

In redundant modeling, each model is logically redundant with respect to the 
other model plus the channeling constraints. In general, the propagators defined 
for two viewpoints act in different ways and discover information at different stages 
in the search. However, we show two possibilities in which propagation caused 
by some constraints in one model can be made redundant by: (a) propagation 
induced from constraints in the other model through channels and (b) propagation 
of the channels themselves. In order to show that, we need some useful lemmas in 
comparing the propagation strength of logically redundant constraints. 

4.1 Useful Lemmas 

It is clear that a constraint C2 that is logically redundant with respect to constraint 
ci is also propagation redundant with respect to ci . 

Lemma 4.1. If\= Dinu Aci ^ C2 then {dsb{ci)} » {dsb{c2)}. 

Proof. Follows immediately from Lemma 2.5. □ 

Typically though a logically redundant constraint C2 is made logically redundant 
by a conjunction of other constraints. It is well known that in general the domain 
(and set bounds) propagation of a conjunction of constraints is not equivalent to 
applying the domain (and set bounds) propagators individually. 

Lemma 4.2. {dsb{ci AC2)} > {dsb{ci), dsb{c2)}. 

Proof. We just show the case for domain propagation of integer variables. Sim- 
ilar arguments apply for set bounds propagation of set variables. Suppose to the 
contrary that {dom(ci), dom{c2)} eliminates a value d from y € {vars{ci Ac2) n V/) 
where d G D{y), but d G do'm{ci A C2){D){y). By definition of propagation solver, 
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there can be no solution 9 which satisfies ci A C2 in £) where 9{y) = d. Hence, 
{dom{c\ A C2)} must eliminate d from y. Contradiction. □ 

But there is a case where propagation of a conjunction is equivalent to propaga- 
tion on the individual conjuncts. 

Lemma 4.3. Let ci and C2 he two constraints sharing at most one integer variable 
X € V/, then {dsb{ci), dsb{c2)} ~ {dsh(ci A C2)}. 

Proof. We have {dsb{ci A C2)} > {dsb{ci), dsb{c2)} by Lemma 4.2. 

To show {dsb{ci), dsb (02)} ^ {dsb(ci A C2)}, suppose 0i is a solution of ci and 
62 is a solution of C2, where ^i(a;) = ^2(3;). Then if 61 G D and 62 € D we have 
that 01 U 62 €z D and (?i U 02 is a solution of ci A C2 . 

Suppose dsb{ci A C2) eliminates a value d from y G Vi- Then there can be no 
solution 6 of ci AC2 in D, where 6{y) = d. Assume w.l.o.g. that y e vars{ci). Sup- 
pose to the contrary that exists solution of ci, 61 G D where 9{y) ~ d. Now if there 
exists a solution 62 & D of C2 where 92{x) — Oi{x) then we have a contradiction. 
Otherwise there is no such 82 hence dsb{c2){D) eliminates the value Oi{x) from the 
domain of x, and hence eliminates 9i as a solution of ci. 

Suppose dsb{ci A C2) adds a value d to inf£)(S') where S* e Vs. Then there 
is no solution 6 of ci /\ C2 in D, where d ^ d{S). Assume w.l.o.g. that 5* € 
vars(ci ) . Suppose to the contrary that exists solution of Ci , Oi £ D where d ^ 0i{S). 
Now if there exists a solution 62 £ D of C2 where 02{x) = 6i{x) then we have a 
contradiction. Hence there is no such 62 thus dsb{c2){D) eliminates the value 6i{x) 
from the domain of x, and hence eliminates Oi as a solution of ci. 

Suppose dsb{ci A C2) eliminates a value d from sup£)(5) where S G Vs- Then 
there is no solution of Ci A C2 in D, where d S d{S). Assume w.l.o.g. that 
S G vars{ci). Suppose to the contrary that exists solution of ci, 9i E D where 
d ^ 6i{S). Now if there exists a solution 62 & D of C2 where 92{x) = 9i{x) then 
we have a contradiction. Hence there is no such 62 dsb{c2){D) eliminates the value 
0\{x) from the domain of x, and hence eliminates Oi as a solution of ci. □ 

Lemma 4.3 also allows us to always state that {dsh{c), dsb{cx)} ~ {dsb{c A Cx)} 
where Cx is any constraint on a single integer variable. 

Note that the result above does not hold when the single variable shared is a set 
variable, because we only apply set bounds propagation. If we did use set domain 
propagators the result readily extends to the case where a single shared variable is 
a set variable. 

Example 4.4. Consider the set constraints ci = S e {{1},{2,3}} and C2 = 
S e {{2}, {1,3}}, then dsb{ci){D) = dsb{ci){D) = D where D{S) = {0..{1, 2, 3}}, 
while dsb{ci A C2){D) is a false domain since ci A C2 is unsatisfiable. Note that 
dom{ci){D) = Di where Di{S) = {{!}, {{2, 3}} and dom{c2)iDi){S) =0. □ 

A corollary of Lemma 4.3 is that we can determine domain consistency of an entire 
integer CSP with tree structure just using the individual domain propagators, since 
we can repeatedly apply the above lemma to break the conjunction of constraints 
into individual constraints. This is quite related to the "backtrack-free" approach 
to solving CSPs with tree structure of Preuder [Preuder 1982]. 
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4.2 Propagation Redundancy Through Channels 

To show that the propagation caused by some constraints in one model is subsumed 
by propagation induced from constraints in the other model through channels, we 
need to break up the consideration of a constraint into individual propagation 
rules. The following lemma ensures that the domain and set bounds propagator of 
a constraint is equivalent to the union of the propagation rules implemented by the 
propagator. 

Lemma 4.5. Given a constraint c. Then {dsb{c)} ~ prop{dsb{c)). 

Proof. We have {dsb{c)} » prop{dsb{c)) by Lemma 2.5 and and Lemma 4.1. 
It remains to show that prop{dsb{cj) ^ {dsb{c)}. Suppose D' = dsb{c){D), and 
D" = solv{prop{dsh{c)),D), We have the following cases: 

— Suppose d ^ dom{c){D){xi) and d G D{xi) where Xi G V/. Then clearly dsb{c) 
implements a rule r = D ^ Xi d where D = AJ^^ /\d'eD,„t{xj)-D{xj) xj ^ d'. 
Now since r G rules{dsb{c)) we have that prop{dsb{c)) {r} by the definition 
of prop{dsb{c)). Hence d ^ D"{xi). 

— Suppose d € info'iSi) and d ^ mf£){Si) where 5, e Vs- Clearly dsb{c) imple- 
ments a rule r = D ^ d G Si where 

Afc=l Ad"6supo_^_^(Sfc)-sup^(Sfc) d ^Sk 

As above since r e rules{dsb{c)) we have that prop{dsb{c)) » {r} by the defini- 
tion of prop{dsb{c)). Hence d G info' (Si). 
— d ^ sup^,(S'i) and d G supjj{Si) where Si e Vs- This case is similar to the 
previous case. 

Since D" has all the modifications to the domain made by D' we have that D" C D' 
and prop{dsb{c)) » {dsb{c)}. □ 

A propagation rule Ci ^ ci directly subsumes a rule C2 ^ C2 if |= {Dimi AC'2) ^ 
Ci and \= {D„iit A ci) ^ C2. A propagation rule that is mapped by a channel 
function to a rule directly subsumed by another propagation rule is propagation 
redundant. We extend channel functions to map conjunctions of atomic constraints 
in the obvious manner •0(ci A • • • A c„) = <^(ci) A • • • A ^(c„). 

Lemma 4.6. Let C ^ a be a propagation rule on Y variables, and C ^ a' be a 
propagation rule on X variables. If 

\= D,mt A 0-^0) ^ C" and h A™t A c' ^ 0'\c), 

then 

{C ^c'}UF^:^{C^c}. 

Proof. Suppose C ^ c is applied to enforce c. Then \= D ^ C. Now if 

Di = solv{F(^, D) then by the definition of clearly \= Di ^ (}~^{C) and hence 
\= Di ^ C . By the definition of the propagation rule C" ^ c' we have |= D2 c', 
where D2 = (C ^ c')(£'i); hence \= D2 ^ <>'^{c)- Then, by the definition of F^, 
1= Ds — > c, where D3 = solv{F^,D2). □ 
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We can straightforwardly lift the above results to talk about propagation rules 
that are directly subsumed by the domain and set bounds propagator for a con- 
straint, and then lift to a set of propagation rules implemented by some propagator. 

Theorem 4.7. Let C c be a propagation rule on Y variables, and cx be a 
constraint on X variables. If 

N (Amt A CX A 0-\C)) ^ 0-\c), 

then 

{rfs6(cx)}UJ'o >{C^c}. 

Proof. By Lemma 4.5 we have that dsb{cx) implements the propagation rule 

<0>~^(C) ^^^(c). Hence by Lemma 4.6 the result holds. □ 

Corollary 4.8. Let cy be a constraint on Y variables, and cx be a constraint 
on X variables. If 

\= {Dinu A A O-i(C)) ^ 0"'(c) for all C ^ c e prop{dsb{cY)), 

then 

{dsb{cx)} U Fo > {dsbicy)}. 

Example 4.9. Consider the (LY2.1) constraints cy = = 3i — 2 <^ = 
3i — 1 of the Langford's Problem (Example 3.1), where 1 < « < 9 and 1 < j < 
27 — 2{i + 1). The propagation rules for dsb{cY) are 

yj = 3i-2 ^ Vj+ii+i) = 3i - 1 

VjHi+i) =32-1 ^ yj =3i-2 

yjy^3i-2 ^ Vj+ii+i) 7^ 3i - 1 

7^ 3i - 1 ^ yj j^3i-2 

We have that for cx = xsi-i = xm-i + (i + 1) of (LX2.1), ^ Dinn A cxA M 
(C) — > 1X1 (c) for all propagation rules above. For example the first rule is mapped 
to x-ii-1 = j ^ xsi-i = j + i + I. Hence the conditions of Corollary 4.8 hold and 
dsb{cY) is propagation redundant. 

Wc can similarly show the (LY2.2) constraints Cy = yj ~ 3i~2 <^ yj+2(i+i) = 3i, 
where 1 < i < 9 and 1 < j < 27 — 2(z + l), are propagation redundant using 
Cx A c^, where c'x = Xsi = Xsi-i + (i + 1) of (LX2.2). Although model Mx 
docs not include a domain propagator for cx A c^, wc can still show propagation 
redundancy since {dsb{cx), dsb{c'x)} « {dsb{cx A c'x)} by Theorem 4.3. Similar 
reasoning applies to show that the (LY3) constraints yj ^ 3i — 2, where 1 < i < 9 
and (28 — 2{i + 1)) < j < 27, are made propagation redundant by cx A c^. □ 

For brevity we shall introduce pseudo atomic constraints x < d equivalent to the 
conjunction x ^ d+ 1, . . . ,x sup^).^.^ {x) and x > d equivalent to the conjunction 
X ^ inf Dir,it (x), . . . ,x ^ d — 1, to discuss the next example. 

Example 4.10. Consider the (BX2) constraints cx = = j)) = 1j of 

the balanced academic curriculum problem (Example 3.5), where 1 < j < n. The 
propagation rules C ^ c for dsb{cx) are 

qj <dAXi^=j A--- AXij=j ^ Xi^ j 
Xii=jA---AXi^=j ^ qj > d 
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for all / = {ii, . . . , id} C {1, . . . , m} and i G {1, . . . , m} — /; and 



qj >dAxi^ ^ j A • • • A Xi^_^ 7^ j ^ a;, = j 
a;i^ ^ j A • • • A a;i^_^ ^ j ^ < d 

for all 7 = {zi, . . . , im-d} C {!,..., m} and i e {!,..., m} — /. Notice that all the 

atomic constraints involving qj arc mapped to themselves by {}, since qj is shared 
by the two models. The rules arc mapped to 

qj < d A ii € Sj A ■ ■ ■ A id € Sj ^ i ^ Sj 

ii & Sj A ■ ■ ■ Aid & Sj ^ Qj < d 

qj > d Ail ^ Sj A ■ ■ ■ A im-d <^ Sj ^ i G Sj 

ii ^ Sj A ■ ■ ■ A im-d Sj ^ Qj > d 

We have that for cs = \Sj\ = q^ of (BS3), [= {D,nu A eg A {}(C)) ^ {}(c) for 
all the propagation rules above. Hence, dsb{cx) is propagation redundant using 
Corollary 4.8. 

Similar reasoning applies to show that each constraint {Y^'^i{{xi = j) x ti)) = Ij 
of (BXl) is made propagation redundant by i^i^g. ti) = Ij of (BS2), where 1 < 
j <n. n 

Often a single constraint does not capture all the propagation effects of a con- 
straint on the other side of the permutation model. In that case we may need to 
find for each particular propagation rule, a constraint on the other side that causes 
the same propagation to occur. 

Theorem 4.11. Let cy he a propagator on, Y variables. Suppose for each r = 
(C ^ c) S prop{dsb{cY)), there exists constraint imp{r) on X variables where 
h (Amf Amp(r)AO-^(C)) 0"^(c), then\Jr(^prop(dsb(cv)){dsh{imp{r))}i^F(^ > 

{dsh{CY)}. 

Proof. The proof follows straightforwardly from Lemma 4.5 and Theorem 4.7. □ 

Example 4.12. Consider the (IY2.1) constraints cy = (jji - jjj = 1) ^ b-'j-i = 
Uj) of the all intervals series problem (Example 3.2), where 1 < i < j < n. The 
propagation rules for dsb{cY) have the forms 

rl : Ui = k + 1 Ayj = k ^ Vj-i = k 
r2 : Vj-i ^ k Ajij = k ^ jji ^ k + 1 
rS : yi = k + lAl^yj^k 

where in r3, / is any conjunction of disequations on Vj-i and yj, not including yj ^ k 
ensuring that Vj^i ^ yj. We can show for imp{rl) = imp{r2) = imp{r3) = {uk = 
\xk - Xk+i\) of (1X2) that 1= (£)„„( A imp{rl) A Xk+i = » A Xfe = j) ^ (uj, = j - i) 
and \= {Dinit A imp{r2) Auk i= j — i A Xk = j) {xk+i 7^ i)- For the remaining 
propagation rules (r3), it is clear that / must contain Vj^i / k since it does not 
contain yj ^ k and it must force the two to be different. We can show that 

h {Dimt A imp{ri) Auk^ j - i A Xk+i = i) ^ {xk ^ j)- 

Note that even though each example propagation rule rl-r3 is made propagation 
redundant by the same constraint Uk = \xk — Xk+i \ we require a diff'erent constraint 
for each diflFerent value of k. 
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Hence the constraint is propagation redundant by Theorem 4.11. Similarly for 
the other (IY2.2) constraints {yj — yj = 1) =4- {vj-i = yi), where 1 < i < j < n. 
Clearly, the rcdimdant constraint (IY3) — y„| = 1) A {vn-i = min(j/i,2/„)) is 
non-propagation redundant. □ 

4.3 Propagation Redundancy Caused by Channels 

The channels themselves may actually restrict the possible solutions in one or both 
models involved. We concentrate on the X model, since the restrictions on the Y 
model can be seen by examining the inverse channel function. 

A channel function <^ is restrictive (on the variables X) if ^ D^nit ^FC^, 
that is not all valuations on X variables are extensible to solutions of Co- 

Example 4.13. The permutation channel function Xl is restrictive, for example 

{xi = 2,X2 = 2} cannot be extended to be a solution of Cm, since it requires 7/2 to 
take both values 1 and 2. The Boolean channel function A is unrestrictive. Any 
valuation on X variables can be extended to a solution of Ca- However is 
restrictive, for example {zn — 1, Z12 — 1} cannot be extended to a solution of Ca 
since it requires xi to be both 1 and 2. Similarly the set channel function {} is 
unrestrictive while {}~^ is restrictive. For example = {1},<S'2 = {1} cannot be 
extended to a solution of C{} since it requires xi to be both 1 and 2. The set2bool 
channel — is unrestrictive in both directions. □ 

4.3.1 Restrictive Channel Functions. Restrictive channel functions can them- 
selves make constraints propagation redundant. 

Smith [2000] first observe that the permutation channel makes each of the dis- 
equations between variables in either model propagation redundant. Walsh [2001] 
proves this holds for other notions of consistency. 

Lemma 4.14 Walsh [2001]. Fx, » {dsb{xi ^ xj)} for alll<i<j <n. □ 

Example 4.15. Lemma 4.14 makes the (LXl) and (LYl) constraints of the 

Langford's Problem (Example 3.1), and the (IXl.l), (1X1.2), (lYl.l) and (IY1.2) 
constraints of the all intervals series (Example 3.2) propagation redundant. □ 

Implicit in the Boolean channel is that each integer variable can take only one, 
and must take one value. This is represented in the Boolean model as the constraint 
Sj=o '^ij ~ ^- enforced by the restrictive channel function A~^. 

Lemma 4.16. Fa > {c?s6(Ej=i = 1)} for alll<i<n. 

Proof. The propagation rules for dsb{J2j=i Zij = 1) are 

= 1 ^ Zij' = 0,jV f 
Zii = 0,..., Zi(^j-i) = 0, =0,...,Zik = ^ Zij = 1 

For the first rule, D{zij) = {1} means that dsb{J2'j=i Zij = 1)(-D)(%') ^ {0}- Let 
D = soIv{Fa,D). Then D{xi) = {j} using Xi = j Zij = 1 and 1 ^ D{zij'),j 7^ f 
using Xi = j' <^ Zij' = 1. 

For the second rule, D{zij) = {0},V1 < j' ^ j < k and dsb{J2j=i Zij = 
l)iD){zij) C {1}. Let D = solv{F^,D). Now D{xi)f){l, . . . + . . . ,k} = 9 

18 



using Xi = j' Ziji = 1. Hence D{xi) = {j} and ^ D{zij) using Xi = j <^ Zij = 
1. □ 

Example 4.17. The (QZl) constraints of the n-Queens Problem (Example 3.3) 
are propagation redundant using Lemma 4.16. □ 

The channel function {} ^ is restrictive, since each variable Xi E X can only take 
a single value j. It means that Sj fi Sy = for all < j < < m. It is clear that 
F{j makes these constraints propagation redundant. 

Lemma 4.18. F^y > {dsb{Sj n Sy = 0)} for all I < j < f < m. 

Proof. The propagation rules for dsb{Sj fl Sj' = 0) are 

% ^ Sj * i ^ Sj' 
i € Sj' ^ i ^ Sj 

For the first rule, assume w.l.o.g. that dsb{Sj Ci Sj' = 0) makes i ^ 5^/, then 
i G infi3(S'j). Clearly, D' = soIv{Fq,D) is such that D'{xi) = {j} using Xi = j ^ 
i G Sj. Then i ^ sup ]j,{Sj), using Xi = j' i & Sj'. A similar arguments applies 
to the second rule. □ 

Example 4.19. Lemma 4.18 makes both the (GSl) constraints of the social 
golfers problem (Example 3.4) and (BSl) constraints of the balanced academic 
curriculum problem (Example 3.5) propagation redundant. □ 

4.3.2 Unrestrictive Channel Functions. Unrestrictive channel functions do not 
make any constraints (on X) propagation redundant. Interestingly in this case we 
can argue about propagation redundancy simply in terms of logical consequence. 

Theorem 4.20. Let {} be an unrestrictive channel function, let cy be a con- 
straint on Y variables, and cx a constraint on X variables. If 

\= {Dinit ACX A C^) CY, 

then 

{dsb{cx)} U F.^ » {dsbicy)}. 

Proof. Let C cG prop{dsb{cY)) then clearly |= D„,it A Cx A Co — »■ (C — > c) , 
since |= Dinu Acy ^ {C ^ c) by Lemma 2.5. We show that 

\= Dir^u ACX ^ {<>-\C) ^ 0-\c)) 

by contradiction. 

Assume to the contrary, there is a solution 9x of Dinu Acx A <C>~^ (C) A -i<C>~^ (c) . 
Now since 6x is a solution of Dinu A Cx and 'C> is non-constraining it is extensible 
to a solution 9 = 9x ^ Oy of . 

Now 9y is a solution of C by the definition of , but similarly since 9x is not 
a solution of <C>~^(c), then 6y is not a solution of c. Contradiction. □ 

The reason the channel function must be unrestrictive for this result to hold is 
that the |= {Dinu Acx A Co) — > cy is too weak a condition in the general case. 
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Example 4.21. The permutation channel function is restrictive. Now \= C ^ 
?/3 = 3, where C = xi +X2 < 4aCm, since the only solutions of C are {a;i i— > 1, X2 i— > 
2,X3 I— > 3,7/1 ^ 1,1/2 ^ 2,1/3 '-^ 3} and {xi 2,0:2 i—> 1,2:3 ^ 3,yi 2,?/2 
1)2/3 3}. But clearly it is not the case that X1+X2 < 4 — > xa = 3. The problem is 
that the channel Cm removes solutions of a;i +a;2 < 4 like {xi 1, a;2 1-^ 1, X3 1— > 1} 
from consideration. □ 

We can use Theorem 4.20 to prove propagation redundancy of many of the prop- 
agators in our examples. 

Example 4.22. Consider the (QZ3.1) constraint cz = J2i=i^ii ^ 1 of the n- 
Queens Problem (Example 3.3). We can show for 

Cx = Xi ^ Xi — i + 1 A ■ ■ ■ A Xi-i ^ Xi — 1 A Xi^i ^Xi + lA---AXn^Xi+n — i 

of (QX2.1) and (QX2.2) that |= Dinit Acx ACa^ cz- Now 

{dsb{cx)} ~ {dsb{xi ^ Xi — i + 1), . . . , dsb{xn Xi + n — i)} 

by Lemma 4.3. Since A is an unrcstrictivc channel function, by Theorem 4.20 we 
have that dsb{cz) is propagation redundant. A similar argument applies to the 
(QZ3.2), (QZ4.1), (QZ4.2), (QZ4.3) and (QZ4.4) constraints. 

Note that the (QZ2) constraints X]"^i Zij = 1, where 1 < j < n, arc not propaga- 
tion redundant, although the constraint Ym=i ^ij < 1 is made redundant by (QXl) 
(using a similar argument to the (QZ3.1) constraints). □ 

Example 4.23. For the social golfers problem (Example 3.4), the (GS2) con- 
straints = \Sik\ = s where 1 < « < .9 and 1 < /c < w, we can show for 
cxi = Sf^^(a;;fc = i) = s of (GXl) that |= Di^t A cx^ A C{} cs^. For 
the (GS3) constraints cs^ = \Si^ki H Si^k2\ ^ 1 where I < ii ^ 12 < 9 and 
1 < ki ^ k2 < w, we can also show for 0x2 ^ ~'(a^/ifei = a^^fci ^ ^-hko = xi.-,k-,) of 
(GX2) that 1= Dinit l^cx-i ^ CS2- Since {} is an unrestrictive channel function 

by Theorem 4.20 we have both dsb{csi) and dsb{cs^) are propagation redundant. 
□ 

Example 4.24. For the balanced academic curriculum problem (Example 3.5), 
the (BS4) constraints cs = {i & Sk) =^ (j ^ Sk') where (i, j) G i?, 1 < fc < n— 1, and 
1 < fc' < fc, we can show for cx = Xi < xj of (BX3) that |= {Dinit Acx AC^y) — * cs- 
Since {} is an unrestrictive channel function by Theorem 4.20 we have that dsb{cs) 
is propagation redundant. □ 

In part because the — channel is unrestrictive in both directions, we can prove 
that set bounds propagation does not provide any more propagation strength than 
the mapping of set constraints to Booleans. First, we introduce the notion of nogood 
constraints. 

A nogood constraint c is one where every valuation over vars{c) in Dinit except 
one valuation ^ is a solution of c. We call the non-solution valuation 9 the nogood 
of c. The following lemma identify the condition where propagation of the conjunc- 
tion of a nogood constraint ci and a constraint C2 of integer variables such that 
vars{c2) C vars{c\) is equivalent to propagation on the individual conjuncts. The 
condition requires that each valuation 6' G Dinit differing from the nogood 9 by 
only one assignment must be a solution of C2. 
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Lemma 4.25. Let ci be a nogood constraint with vars{ci) = {xi, . . . ,a;„} C V/ 
and the nogood 6, and ci he a constraint with vars{c2) C vars{c-i). If for all 

valuations 6' G Dinit, such that there exists 1 < j < n and 9'{xi) = 9{xi) for all 
1 £ * 7^ J O'TC solutions of C2, then {dom{ci), dom{c2)} ~ {dom{ci A C2)}. 

Proof. By Lemma 4.2 we have that {dom{ci A C2)} ^ {dom{ci), dom{c2)}- 

It remains to show {dom(ci), dom{c2)} {rfom(ci A C2)}. Suppose to the con- 
trary that d ^ dom{ci /\C2){D){xk) and d € solv{{dom{ci), dom{c2)}, D){xk) where 
1 <k <n and d G D(xk)- 

First, we show that the nogood 6 G D. Otherwise, when 6 ^ D, there must exists 
1 < J < such that 0{xj) ^ D{xj). Then any vahiation 0" G D on vars{ci) which 
makes C2 true is also a solution of Ci since 9"{xj) / 6{xj). Thus dom{c2) removes 
any values removed by dom{ci A C2), Contradiction. 

Second, we show that 9{xk) = d. Otherwise, when d{xk) 7^ d, any solution 9" G D 
on vars{ci) which makes C2 true is also a solution of ci. Hence d ^ dom{c2){D){xk). 
Contradiction. 

Now, since d G dom{ci){D){xk) , it is not the case that D{xi) = {9{xi)}, for 
1 < k ^ i < n. Thus there exists j ^ k and Xj G wars(c2) such that \D{xj)\ > 2. 
Then we have dj G D{xj) such that 7^ 9{xj). 

Consider the valuation 9' defined as 6''(xi) = 9{xi),\ < i ^ j < n and 9'{xj) = dj, 
note that 9'(xk) = 9(xk) = d. By construction 9' € D is a solution of ci and also 
of 02 by the conditions of the lemma. Hence, d G dom{c\ A C2)(-D)(a;fe) by the 
definition of dom. Contradiction. □ 

The following theorem prove that set bounds propagation of set constraints is 
equivalent to domain propagation of the corresponding Boolean constraints. It does 
however still provide a more efficient implementation. 

Theorem 4.26. Let dsb{c) be the set bounds propagator for set constraint c, let 
^c) be the Boolean equivalent of c, then (a) {dsb{c)}U ^ {dsb{c') \ c' G —(c)} 
and (b) {dsb(c') \ d G -(c)} UF^ > {rfs6(c)}. 

Proof. Wc can verify that |= D„att\C^ ^ (c ^ [l\{c! \ c' G —(c)})) for each of 
the constraints c above. Since — is an unrcstrictive channel this immediately gives 
us that {dsb{c)} \JF^ ^ {dsb{c') \ c' G —(c)} using Theorem 4.20 and Lemma 4.1. 

Now since — is also an unrcstrictive channel we also have by Theorem 4.20 that 
{dsb{^{c' I c' G -(c)})} UF^ > {dsb{c)}. It remains to show that {dsb{h{c! \ d G 
-(c)})} « {dsb{c') I c' G -(c)}. 

For c of the form Si — %, Sa ^ Sb^ SaC\ Sb = %, and |5'i| = m no two constraints 
in —(c) share a variable. Hence the result holds by Lemma 4.3. 

For the remaining constraints wc can, again by Lemma 4.3 simply consider the 3 
Boolean constraints for a particular j, since the remaining constraints do not share 
variables. Wc show the case for c = Sa = SbUSc, the other constraints Sa = Sbr\Sc 
and Sa = Sb — Sc are similar. Wc first show that 

{dsb{Zaj < Zbj + Zcj A {Zbj < Zaj A Zcj < Zaj)} 
W {dsb{Zaj < Zbj + Zcj), dsb{zbj < Zaj A Zcj < Zaj)} 

Now Ci = Zaj < Zbj + Zcj is a nogood constraint with nogood {zaj ^,Zbj 1— *■ 
0,Zcj I— > 0}, while the valuations {zaj 1— > 0,Zbj 1— > 0,Zcj 0}, {zaj l,^6j * 
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1, Zcj I— 0}, {zaj ^ 1, Zbj ^ 0, Zcj ^ 1} are all solutions of C2 A C3 where C2 = < 
Zaj and C3 = A^cj < -Zaj- Hence by Lemma 4.25 we have that {dsh{c\ A (c2 AC3))} « 
{rfs6(ci), dsh{c2 A C3)}. 

Now C2 and C3 share only one variable, hence {dsh{c2 A C3)} ~ {ds6(c2), ds6(c3)} 
by Lemma 4.3. Hence {dsh{c\ A (c2 A C3))} » {rfs6(ci), dsb{c2), dsb{cs)} as de- 
sired. □ 

5. DISCUSSIONS 

Propagation redundancy reasoning is highly sensitive to the exact form of the con- 
straints in the model. A slightly different model may have different degree of 
propagation redundancy. 

Example 5.1. Another version of the second model of the all intervals series 
constraints, replaces (IY2.1) (y^ — yj = 1) {vj-i = yj) and (IY2.2) {yj — = 

1) =^ (^j-j = Vi) by the (IY4) constraints 

QVi - = 1) ^ {vj-i = min(t/i,i/j)), 

where 1 < i < j < n. This combined constraint is not propagation redundant. 
There are cases where it can propagate information not discovered by the first 
model. For example cy = {\y3 — yr] = 1) => {v4 = m.m{y3,yr)) when D{y3) = {5} 
and D{yr) = {4,6}, makes ds6(cy)(D)(w4) = {4,5}, but there is no equivalent 
propagation in the X model. □ 

Our definition of set bounds propagators agree with the definitions of the original 
set bounds propagation system [Gervet 1997]. There are stronger set based prop- 
agators (see [Azevedo and Barahona 2000; Miiller 2001]) that reason more about 

cardinalities. 

Example 5.2. Suppose we have constraints Si C ^2, = 2, |S'2| = 1 where 
Si and 52 range over {1,2,3}. The propagators {dsb{Si C S2)} U {ds6(|S'i| = 

2) } U {dsb{\S2\ = 1)} lead to no propagation. But if we propagate cardinality 
information, wo immediately get failure. □ 

We can model most common cardinality reasoning using additional propagators. 
Let every set variable S has an attached integer variable xs representing its cardi- 
nality. This attachment comes with implicit cardinality propagator: dsh{\S\ = xs)- 
We can then understand stronger cardinality reasoning as additional constraints 
using these implicit cardinality variables. 

Example 5.3. With cardinality reasoning, the set constraints c = S C S' de- 
fines implicitly the constraint: xs < xs' where \S\ = xs and \S'\ = xs'- Hence, the 
correct propagators for c are {dsb{S C S")} U {dsb{xs < xs')}, together with the 
implicit propagators, {dsb{\S\ = xs)} U {dsb{\S'\ = xs')}- □ 

With cardinality reasoning, for a set constraint to be propagation redundant, 
we need to prove the redundancy of the extra cardinality propagators too. For the 
examples in this paper, this is straightforward. However, some cardinality reasoning 
is difficult to model using additional propagators as the reasoning relies on dynamic 
information during propagation. 
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Example 5.4. The set constraint x = J2ies'^i where U = and 
Dinit{S) = {0 . . . {1, . . . , n}} is supported in ILOG Solver 4.4. The cardinality 
reasoning can be computed efficiently by: 

x>{ Yl + (mm{«,|i e sup(5)} * {\S\ - I inf(5)|)) 

ieinfc(S) 

and 

a; < ( V Ui) + {ma^{uj\j e sup(5)} * - | inf(5)|)), 
— ^ r> ^ 

where D is the current domain during propagation. Such reasoning is rather difficult 
to model as separate propagators. □ 

6. EXPERIMENTS 

We can take advantage of the reasoning about propagation redundancy to eliminate 
propagators that arc propagation redundant. We then get a model with exactly 
the same propagation strength but with less propagators. This can translate into 
faster propagation.^ 

In the following experiments, all the benchmarks are executed using ILOG Solver 
4.4 on Sun Ultra 5/400 workstations running Solaris 8. The first column of each 
table describes the models under comparison. In the case of combined models, we 
have the choice of labeling the variables from just one model, or from both models 
together. The second cohimn indicates the choices of search variables. However, 
the question of choosing the "best" set of search variables that gives the smallest 
search space is out of the scope of this paper. The third column and beyond give the 
problem instances. For each problem instanc;e, we measure the total number of fails 
and CPU time (in seconds) to compare the performance of the different models. 
In other words, each row shows the performance of a particular model across the 
different problem instances. Table entries marked with a " " mean failure to solve 
the problem after one hour of execution. Table entries in bold letters highlight the 
model with the smallest number of fails and best runtime for each problem instance. 

6.1 Langford's Problem 

Table I compares the different models for finding all the solutions of the Langford's 
Problem. Variable ordering is based on the first-fail principle and value ordering 
is from the least to the greatest. The models under comparison include the single 
models: Mx and My, the full combined model Mx + Cm + My, and an optimized 
combined model LX2.1 + LX2.2 + Cm as discussed in Examples 4.9 and 4.15. The 
opt model corresponds to the minimal combined model of Smith [2000] . Our results 
agree with those presented by Smith, where the opt models are faster and maintain 
the same number of fails as the full model for all three set of search variables. The 
opt model with search variables X UY is the fastest among all the models under 
comparison. 



*Note there is no guarantee since e.g. the number of propagation steps may have increased. 
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Table I. Results of the Langford's Problem 





Seareh 


(3 X 10) 


(3 X 11) 


(4 


X 14) 


(4 X 15) 


Model 


Variables 


fails (sec) 


fails (sec) 


fails 


(sec) 


fails (sec) 


Mx 


X 


3114 (1.33) 


14512 (6.09) 


83068 


(87.21) 


351126 (430.54) 


full 


X 


1318 (13.66) 


5177 (62.5) 


20885 (1499.09) 


78556 (6264.74) 


opt 


X 


1318 (1.21) 


5177 (4.77) 


20885 


(51.1) 


78556 (183.38) 


My 


Y 












full 


Y 


1059 (11.17) 


3958 (50.67) 


8139 


(707.77) 


25270 (2450.41) 


opt 


Y 


1059 (0.94) 


3958 (3.49) 


8139 


(22.33) 


25270 (69.19) 


full 


XVJY 


768 (9.56) 


2952 (42.83) 


6553 


(643.04) 


20526 (2293.2) 


opt 




768 (0.77) 


2952 (3.08) 


6553 


(17.62) 


20526 (57.52) 



6.2 All Interval Series 

Tabic II compares the different models for finding all the solutions of the all inter- 
val series problem. Variable ordering is based on the first-fail principle and value 
ordering is from the least to the greatest. The models under comparison include 
the single models: Mx and My, the full combined model Mx + Cm + My, and 
an opiimized combined model 1X2 -|- Cm -|- IY3 as discussed in Examples 4.12 
and 4.15. Puget and Regin [2001] show that all the solutions can be found more 
efficiently by replacing (IXl.l) and (1X1.2) by (IXl.l') and (1X1.2') respectively 
where, rather than using a set of separate disequality constraints, we use just two 
alldif f erent global constraints 

—(IXl.l') alldif ferent([a;i, . . . , a;„]) 
—(1X1.2') alldif ferent([ui, . . . , u„-i]) 

The propagator ds6(alldif f erent([ui, . . . , Vk])) is equivalent to 

ds6(Ati'A^t,+i«»7^«,) 

and has an effective implementation [Regin 1994]. The pr model uses IXl.l' + 
1X1.2' + 1X2. The pr full model is the combination of pr and My, IXl.l' + 1X1.2' 
-I- 1X2 -I- Cm + My- The pr opt model is the optimized combination of pr and 
My, IXl.l' + 1X1.2' + 1X2 + Cm + IY3 since by Lemma 4.1 the same redundancy 
reasoning applies. 

Compared with the single models Mx and My, clearly the full and pr full models 
reduce the number of fails significantly (given the same set of search variables as 
the single models). The opt and pr opt models maintain the same number of fails 
as the full and pr full models respectively for all three set of search variables. Using 
search variables Y, the opt model is the fastest for the smaller instance 12 and the 
pr full model is the fastest for larger instances 13, 14 and 15, as the alldif f erent 
constraints is computationally too expensive for the smaller instance. Note that 
the optimized models opt and pr opt with search variables Y can solve the size 15 
instance much faster than pr, and no other models can solve this instance within 
the time limit. 
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Tal>lc 11. Results of Ihc All liilcrval Series l'rol>lcm 



Model 


Search 
Vars 


n = 12 
fails (sec) 


n= 13 
fails (sec) 


n= 14 
fails (sec) 


n= 15 
fails (sec) 


pr 


X 


38778 (24.32) 


156251 (105.26) 


674346 (530.47) 


3045037 (2328.57) 


Mx 
full 
opt 
pr full 
pr opt 


X 
X 
X 
X 
X 


880112 (260.92) 
39241 (222.07) 
39241 (36.34) 
38461 (236.42) 
38461 (42.77) 


4914499 (1589.83) 
158368 (1048.19) 
158368 (157.84) 
155183 (1088.91) 
155183 (188.94) 


685301 (770.57) 
670045 (910.90) 




My 
full 
opt 
pr full 
pr opt 


Y 
Y 
Y 
Y 
Y 


16280 (70.81) 
16280 (6.36) 
12296 (62.96) 
12296 (7.91) 


62949 (303.61) 
62949 (26.00) 
43681 (260.90) 
43681 (25.78) 


266130 (1458.74) 
266130 (108.54) 
164841 (1127.64) 
164841 (101.42) 


1275661 (553.45) 
704097 (458.12) 


full 

opt 
pr full 
pr opt 


XUY 
XUY 
XUY 
XUY 


39195 (222.42) 
39195 (36.36) 
38447 (230.65) 
38447 (42.47) 


158282 (1065.77) 
158282 (158.40) 
155176 (1094.61) 
155176 (198.36) 


684592 (783.01) 
669950 (898.66) 





6.3 n-Queens Problem 

Table III compares the different models for finding all the solutions of the n-Queens 
problem. Variable ordering is based on the first-fail principle and value ordering 
is from the least to the greatest. The models under comparison include the single 
models: Mx and Mz, the full combined model Mx + Cix + Mz, and an optimized 
combined model QXl + QX2.1 + QX2.2 + Ca + QZ2 as discussed in Examples 4.17 
and 4.22. In fact we can show that only the part ^ij > 1 of the (QZ2) 

constraints in the opt model is not propagation redundant, and we refer to this 
further optimized model as the part model. 

Although the full model with search variables X docs reduce the mimbcr of fails 
with respect to the Mx model, but the expensive propagation of the Mz model 
overwhelms the benefits. For the other two sets of search variables, the full model 
does not reduce the search space. The opt models have the number of fails as the full 
models for all three set of search variables, but we only gain slight improvement 
in runtime for search variable X, we get worse runtime with search variable Z 
and mixed results with search variable X U Z. This is because we require more 
propagation steps to obtain information on the Z variables which are driving the 
search in both cases. Compared the part models against the opt models for all 
three set of search variables, the number of fails are still the same, but the results 
for runtime are mixed. This is again due to change in propagation steps after the 
Boolean model is altered. The single model Mx is still the fastest to find all the 
solutions, followed by the single model Mz- The fact that combining integer model 
with Boolean model incurred too high runtime overheads, even after removing the 
redundant propagators, make the combined models impractical. Note that the n- 
queen problem can be solved more efficiently by combining AIx and My using the 
permutation channels^ together with a clever value ordering heuristic [Cheng et al. 



^The My model of the n-Queens problem has the same set of constraints as Mx- 
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Table III. Results of the n-Queens Problem 



Model 


Seareh 
Variables 


n=ll 
fails (sec) 


n = 12 
fails (sec) 


n = 13 
fails (see) 


n = 14 

fails (sec) 


Mx 
full 
opt 
part 


X 
X 
X 
X 


21796 (2.30) 
17601 (9.43) 
17601 (8.46) 
17601 (8.46) 


101882 (11.27) 
80011 (45.37) 
80011 (42.34) 
80011 (42.58) 


515238 (59.55) 
392128 (240.28) 
392128 (225.30) 
392128 (236.02) 


2830370 (324.95) 
2101047 (1350.05) 
2101047 (1267.51) 
2 10 1047 (1273.12) 


Mz 
full 
opt 
part 


z 
z 
z 
z 


23515 (3.54) 
23515 (9.63) 
23515 (10.40) 
23515 (10.48) 


111076 (15.18) 
111076 (50.04) 
111076 (53.74) 
111076 (57.33) 


561362 (86.12) 

561362 (274.43) 
561362 (307.85) 
561362 (295.14) 


3079792 (482.95) 
3079792 (1591.31) 
3079792 (1715.47) 
3079792 (1706.97) 


full 
opt 
part 


xyjz 
xuz 
xuz 


19609 (9.90) 
19609 (9.93) 
19609 (9.56) 


90933 (50.89) 
90933 (48.69) 
90933 (49.12) 


453886 (274.82) 
453886 (276.30) 
453886 (266.91) 


2463621 (1837.08) 
2463621 (1521.17) 
2463621 (1685.25) 



Table IV. Results of Social Golfers Problem 



Model 


Search 
Variables 


4-3-4 
fails (sec) 


7-2-13 
fails (see) 


8-4-9 
fails (sec) 


9-2-17 
fails (sec) 


Mx 

full 
opL 


X 
X 
X 


1509146 (374.29) 
1509116 (818.20) 
1509110 (758. ()1) 


63860 (34.51) 
63860 (110.83) 
():-!8()0 (51.81) 


32 (0.37) 

32 (1.01) 
32 (0.50) 


7355 (7.89) 

7355 (33.01) 
7355 (11.55) 


Ms 
full 
opt 


S 
S 
S 


2389 (0.43) 
1102 (0.56) 
1102 (0.41) 


37158 (34.91) 
27998 (42.45) 
27998 (17.74) 




74098 (130.72) 
51444 (145.94) 
51444 (57.97) 



1999]. 

6.4 Social Golfers Problem 

Table IV compares the diflFerent models for finding the first solution to the social 
golfers problem. For this problem, we refrain from finding all solutions due to the 
large number of symmetric solutions. Variable ordering for variables X is based on 
players by weeks, and variable ordering for variable S is based on weeks by groups. 
Value ordering is from the least to the greatest for both X and S. The models under 
comparison include the single models: Mx and Ms, the full combined model Mx 
+ C{} -I- Ms, and an optimized combined model GXl + GX2 + C^y as discussed 
in Examples 4.19 and 4.23. The problem instances use the parameters g-s-w as 
described in Example 3.4. 

The full model with search variables S reduces the number of fails when compare 
to Ms, but not the case when compare to Mx with search variables X. The opt 
models speed up the search while keeping the same number of fails as the full 
models for both set of search variables. Using search variables X, the search does 
not benefit from the additional propagation of the combined models full and opt. 
In terms of runtime, no one model dominates the others. The opt model with search 
variables S is the fastest for the instances 4-3-4 and 7-2-13, while the single model 
Mx is the fastest for instances 8-4-9 and 9-2-17. 

26 



Table \' 


. Kcsulls 




cadciiiie C'liir 


culuui i^roblcin 




Search 


8 Periods 


10 i'eriods 


12 Periods 


Model 


Variables 


fails (sec) 


fails (sec) 


fails (sec) 


CPLEX 


n/a 


n/a (1.80) 


n/a (2.27) 


n/a (20.32) 


hybrid 


X 


101 (0.61) 


468 (2.20) 


58442 (146.47) 


hybrid 


Boolean 


219 (0.76) 


277 (1.03) 


315 (2.09) 


Mx 


X 


101 (0.04) 


468 (0.25) 


33602 (11.62) 


full 


X 


101 (0.24) 


470 (1.80) 


33530 (192.62) 


opt 


X 


101 (0.08) 


470 (0.68) 


33530 (38.54) 


Ms 


s 








full 


s 


1.577 (2.8.'?) 


32.3 (0.81) 


882 (4.56) 


opt 


s 


1577 (0.94) 


323 (0.24) 


882 (0.95) 



6.5 Balanced Academic Curriculum Problem 

Table V compares the model for finding the optimal solution and proving optimality 
for the problem instances of the balanced academic curriculum problem posted in 
CSPLib. Variable ordering for X is based on the first-fail principle, and variable 
ordering for S is based on ascending order of the indices. Value ordering is from the 
least to the greatest for both set of variables. The models under comparison include 
the single models: Mx and Ms, the full combined model Mx + C^} + Mg, and 
an optimized combined model Bl.l + B1.2 + B2.1 + B2.2 + BX3 + C{} + BS2 + 
BS3 as discussed in Examples 4.10, 4.19 and 4.24. Hnieh et al. [2002] report that 
it is diSicult to find the optimal solution and prove optimality with propagation- 
based constraint solver alone. However, by adding redundant constraints (B2.1) and 
(B2.2), we were able to solve all the; problem instances with Mx alone. The row 
CPLEX gives the runtime for solving the problem instances with ILOG CPLEX 8.0 
using an integer linear programming (ILP) model by Hnich et al. [2002]. The row 
Hybrid implements the hybrid ILP and CP model described by Hnieh et al. [2002] 
together with the redundant constraints (B2.1) and (B2.2) using ILOG Hybrid 1.3. 

Clearly, the full model with search variables S substantially reduces the search 
space when compared with the single model Ms, but mixed results when compared 
to the single model Mx with search variables X. In solving optimization problems 
with propagation-based constraint solver, constraints for reducing the bounds of 
the objective function are added dynamically to reduce the search space. Hence, 
additional propagation may change the ordering sequence of the solutions, which 
can increase or decrease the numbers of fails. This explains why the full model 
with search variables X has more fails than the single model Mx for instances of 
10 Periods. The opt models are faster than the full models while keeping the same 
number of fails for both set of search variables. In terms of runtime, no one model 
dominates the others. The opt model with search variables S is the fastest for 
solving the instances with 10 and 12 periods, while the Mx model is the fastest 
for solving the instances with 8 periods. Although the hybrid models give the least 
number of fails, it suffers from the overhead of invoking two solvers. 
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7. RELATED WORK 



Smith [2000; 2001] has examined the redundant models for a number of individual 
problems including the n-Queens problem, Langford's problem and the social golfers 

problems. She empirically demonstrates that some constraints in the rcdimdant 
models can be removed without increasing the search space. She points out that 
for these problems the so-called minimal combined model, which combine the first 
model and only the variables of the second model (without the constraints) using 
channeling constraints, produces the same search behavior as combining the models 
in full. This is proved in an ad hoc manner by Choi and Lee [2002]. In this paper, 
we aim for a theoretical framework which can determine propagation redundancy 
of a particular constraint involved in redundant models a priori. 

Apt and Monfroy [2001] develop "membership rules" as a way of building propa- 
gators for any constraints. Propagation rules are similar to the "membership rules" 
when restricted to integer variables. However, we develop propagation rules as a 
method for reasoning about the parts of a propagator's behavior. 

Brand [2003] gives a general theorem to determine when a rule is propagation 
redundant with respect to a set of rules in rule-based constraint programming, and 
illustrates the applicability using "membership rules." In fact, our definition of a 
propagation rule satisfies the required properties of Brand's theorem. Hence, we 
can apply Brand's theorem to determine when a propagation rule is propagation 
redundant with respect to a set of propagation rules. In this paper, we are in- 
terested in propagation redundancy beyond the individual propagation rules, but 
propagation redundancy of the constraint as a whole. We also generalize the notion 
of propagation redundancy of a propagation rules through a channel function. 

Hnich et al. [2004] and Walsh [2001] introduce the notion of constraint tight- 
ness as a measure to compare the propagation strength of different permutation 
constraints. Their work focuses on comparing the propagation strength of the 
different notions of consistency over the disequations, channeling constraints, and 
alldif f erent constraints in redundant modeling of only permutation problems 
and injection problems. Our comparison measure is similar to constraint tightness 
except that constraint tightness is parameterized by a local consistency property. 
However, in existing constraint solvers, there are propagators which implement none 
of the (established) local consistency properties. An example is the multiplication 
constraint x = yx z over integer domain as discussed in Apt [2003, pages 219-220]. 
In such cases where the local consistency property of a constraint is unknown, our 
comparison measure woiild still be applicable. In this paper, we are not only in- 
terested in studying the propagation of the permutation constraints, but also the 
other constraints in redundant models. We also cover a broader class of channeling 
constraints beyond the permutation channels. 

Walsh [2003] proves that "bounds consistency" on set (multiset) variables is 
equivalent to bounds consistency on the corresponding occurrence representation. 
This result is related to Theorem 4.26 since the occurrence representation of set 
variables corresponds to Boolean variables described in Section 3.4. However, ex- 
isting constraint solvers break Boolean constraints into parts and propagate each 
part separately. We prove the theorem based on such realistic assumption. 
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8. CONCLUSION 



The contributions of this paper are three-fold. First, we define channehng con- 
straints in terms of channel functions which allow us to cover a broad form of 
redundant modeling. By breaking up a propagator into individual propagation 
rules, we reason that constraints in one model can be made propagation redun- 
dant by constraints in the other model through channels. Second, we introduce the 
notion of restrictive and unrestrictive channel functions to characterize channel- 
ing constraints. Restrictive channel functions can themselves make a constraint in 
the combined model propagation redundant. Unrestrictive channel functions allow 
the detection of propagation redundancy of a constraint in one model with respect 
to a constraint in the other model plus the channels simply in terms of logical 
consequence. Third, benchmarking results confirm that removals of propagation 
redundant constraints from combined model can lead to a faster implementation 
with the same search space. As explained in Section 7, this paper extends related 
work by covering a broader form of redundant modeling and reasoning about the 
propagation redundancy of all the constraints in the redundant models. 

Our work prompts a number of important future directions for research. It is 
interesting to investigate if the process of removing propagation rcchmdant con- 
straints can be (semi-)automated. To use Theorem 4.11 we can straightforwardly 
define the propagation rules for many constraints (parametrically in Dinit ) or con- 
struct them automatically using the approach of Abdennadher and Rigotti [2002] . 
The number of propagation rules for most constraints, however, are exponential. A 
naive implementation is computationally impractical beyond textbook applications. 
A possible approach is to consider parameterized propagation rules, which denotes 
a set of propagation rules, so that the number of rules is vastly reduced. We can 
also try to use Theorem 4.20 to prove propagation redundancy without considering 
propagation rules. 

Variable and value ordering heuristics can affect the amount of search space re- 
duction caused by constraint propagation. There are cases shown in Section 6 in 
which stronger propagation does not reduce the search space. The interesting ques- 
tion is to study a notion of search redundancy which allows even non-propagation 
redundant constraints to be removed without increasing the search space given 
specific variable and value ordering heuristics. 

Redundant modeling gives rise to the need to decide which variables to label dur- 
ing search. As demonstrated in Section 6, the choice of search variables can greatly 
affect the size of the search space. Many related work, for example Geelen [1992], 
Smith [2000; 2001] and Hnich et al. [2004], also shown that certain choices of search 
variables do lead to a smaller search space. Therefore, it is interesting to study and 
establish criteria in choosing the better set of search variables. 
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